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The effect of preheating on the ablative Rayleigh-Taylor instability 
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Abstract. - The evolution of the ablative Rayleigh-Taylor instability (ARTI) in the presence of 
preheating is investigated. When the linear ARTI growth rate has a sharp maximum, a secondary 
spike-bubble structure at a harmonic of the fundamental mode will be excited. Interaction of this 
secondary structure with the latter will eventually lead to its rupture. On the other hand, for 
stronger preheating the linear growth rate has no sharp maximum and in the nonlinear evolution 
no secondary spike-bubble is generated. Instead, the fundamental spike-bubble structure evolves 
into an elongated jet. These results are in agreement with that from the numerical simulations 
and observations. 



During intense-laser interaction with a target in iner- 
tial confinement fusion, a subsonic heat-conduction wave 
is excited at the laser-absorption spot on the target sur- 
face. The intense heat wave propagates into the cold 
dense plasma and generates a shock wave. As the lat- 
ter propagates, it compresses, heats, and accelerates the 
local plasma. A rarefaction wave that decelerates the ac- 
celerating fluid is formed ahead of the heat front, where 
the Rayleigh-Taylor instability (RTI) can take place 1 3 
The RTI is also common in gas combustion, the implosion 
phase of laser-target interaction, as well as in astrophys- 
ical gas clouds [THE]. Linear theory shows that ablation 
caused by the heat front can stabilize the RTI as well as 
the subsequent Kelvin-Hclmholtz instability [ThTTJ. 

Most earlier studies of the ablative RTI (ARTI) used the 
fluid approach and invoked the Spitzer-Harm (SH) con- 
ductivity for the thermal transport. On the other hand, re- 
sults from Fokker-Planck (FP) simulations showed [T^IRTS"] 
that under ICF conditions the SH theory is inadequate 
because of intense plasma heating by the laser-generated 
hot electrons, heat and shock waves, etc. 012], and the 
preheated region in front of the temperature gradient of 
the heat-wave front is often spatially extended. The heat 
transport process involved in the region is rather com- 
plex and several models have been introduced 



Simulations and analytical modeling have shown that the 
preheating can significantly affect the linear growth rate 
and cutoff wavelength of the ARTI, as well as its nonlinear 
development [12, 15, 16 . However, although the onset and 
initial evolution of ARTI have been widely investigated 
and well understood, a more detailed physical understand- 
ing of the later stages of the development is still lacking 

m- 

Here we consider the evolution of the ARTI accounting 
for the preheating effect by analytical modeling and nu- 
merical simulation. The conservation equations including 
laser absorption and electron conduction effects are solved 
numerically using the LARED-S code [17]. The energy 
equation is given by [TH91IT7] 



d t E + V • [(E + P)u] = -V ■ q(T) + S L , 



(1) 



where p, u, T, P and E are the fluid density, velocity, 
temperature, pressure, and total energy E = p(e + u 2 ), 
respectively, e is the specific internal energy, Sl is the 
energy source due to absorption of laser energy by in- 
verse bremsstrahlung, and q = -kVT is the heat flow. 
To include preheating and heat flux limitation effects, we 
invoke the modified SH electron heat conductivity [7 9 
k = nsnf(T), where ksh is the classical SH conductivity, 
and f{T) = 1 + aT" 1 + /3T~ 3 / 2 is the preheat function 
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[5]. Here, T is in units of MK, and a and j3 are free 
parameters of order unity and can be determined empiri- 
cally by comparing with experiments, theories, and/or nu- 
merical simulations. Thus, for constant pressure one can 
define a characteristic preheat length L — T(d x T)^ 1 ~ 
p{dxP)~ l = L SH f(T), where L SH = Ksal PaV a Cp is the 
classical SH length and p a , v a , c p are the density, speed, 
and specific heat, respectively, of the plasma in the abla- 
tion front region. Other heat conduction models [13,18 
can also be used, provided that preheating effects and heat 
flux limitation are properly included. 

In the simulation, we first obtain a quasistationary state 
by irradiating a 200/im thick CH foil target with a 0.35/im 
laser beam at intensity / = Iq sin 4 (t/5ns)/sin 4 (1.0) for 
t < 5 ns and I = I = 10 15 W/cm 2 for t > 5 ns. Such a self- 
consistent one-dimensional quasistationary state appears 
at t > 7 ns, as shown in the left panel of fig. [1] for the 
density, temperature, and pressure near the ablation front. 
To initiate the ARTI on the ablation-front interface, we 
start at this instant a small perturbation given by the 
fluid displacement x = £(y, 0) = 5£ocos(ky), where <5£o is 
the amplitude and A = 2ir/k is the wavelength. In both 
the dense pristine plasma and the ablation regions, 1200 
fine meshes of widths of 0.06 - 0.07 /im are used in the 
x and y directions. The mesh width in the x direction 
is gradually increased on both sides of the region. Up to 
1400 x 200 meshes are used in the simulation. 

The left panel of fig. Q] shows that near the ablation front 
the temperature and density are ~ 0.2 MK and 5.6 - 6.5 
g/cm 3 , respectively, and in the ablation layer they are 
~ 2 MK and ~ 0.8 g/cm 3 , respectively. The right panel 
of fig. Q] shows that the linear ARTI growth rate depends 
rather strongly on the preheating length L. It is thus con- 
venient to define weak preheating (WP) by the condition 
L/Lsn = f(T) <C /o ~ 20, as is the case represented by 
the black squared curve (which is for f(T) ~ 8 at T = 0.2 
MK). Here the linear growth rate has a relatively sharp 
maximum at A ~ 5/im. We define strong preheating (SP) 
by the condition L/Lsn > /o ~ 20. It corresponds to the 
regime where the linear growth rate does not have a sharp 
maximum, or the red circled and green triangled curves. 

We first consider the WP regime. For a — 0.7 and (3 = 
0.4 the preheating length L is about 0.58 /im, with / « 8. 
The initial perturbation has (5£o=0.4/im and A = 12/im. 
Figure Q] shows that the linear growth rate has a sharp 
maximum at X m — 0.375Ai and a cutoff wavelength at 
A c = 0.14Ai. As expected, the linear growth rate agrees 
well with that from the modified Lindl-Takabe formula [7] 
7 = [Akg/0- + a fci)] 1/2 - Pov a k, where A = 0(1) is 
the Atwood number, g is the fluid acceleration, and the 
free parameters cto and j3q are in the range 2-3. That 
is, in the short wavelength regime the effect of ablation 
is significant and the linear growth rate decreases as the 
preheating length L increases. 

The simulation results show that besides the fundamen- 
tal mode (Ai = 12//m), at least six harmonics having wave- 
lengths larger than A c = 0.14Ai appear. The evolution of 
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Fig. 1: (Color online.) Left panel: profiles of the density (0.8 — 
10.5 g/cm 3 ), temperature (0.1 - 2 MK), and pressure (~ 70 
Mb) for the WP regime at t = 9.09 ns near the ablation front 
at x ~ —536. 5/im. The laser is incident from the right. Here 
electron heat conduction heats the overdense region, and the 
ARTI interface is near the density peak. Right panel: linear 
(i.e., initial) growth rate of the ARTI versus the wavelength, 
from numerical simulations. The black curve with squares is 
for the WP case, with a — 0.7 and j3 — 0.4. The red curve 
with circles is for SP1 (a = 5.0, /3 = 1.0), where A c ~ 3/mi, 
L ~ 1.35pm, v a ~ 0.80/im/ns. The green curve with triangles 
is for SP2 (a = 7.0, /3 = 1.4.) 



the fundamental spike £(y — Ai/2, t) is shown in fig. [5] and 
the evolution of the harmonic modes £«(£) at y = Ai/2 for 
n = 1 — 4 are shown in fig. [3] At t = 7.78 ns, the fun- 
damental spike has evolved into a wide spatial structure 
with a bubble of size ~ 0.24Ai. In the plasma region that 
is pre-heated by the heat conduction wave, it grows at 
> 10 6 cm/sec. In this region the temperature (> 1 MK) 
and pressure (> 70 Mb) are also high, but the density de- 
creases sharply, as can be seen in the left panel of fig. [TJ 
The front of the fundamental spike becomes flattened, but 
no mushroom structure is generated. This behavior can 
be attributed to the fact that there is enhanced energy loss 
due to increased electron heat conduction, so that velocity 
shear (and thus fluid rotation) is suppressed [3 . 

On can see that at t = 8.1 ns the fundamental spike 
is still extending forward, but now at the slower rate 
~ 2.5 x 10 6 cm/sec. The displacement £ of the layer is 
close to the modulation wavelength Ai. That is, higher 
harmonics and mode coupling effects become important, 
and the ARTI can be considered to enter a new nonlin- 
ear phase. Here the boundary region of the high-density 
area is gradually ablated because of the high temperature 
there. In fact, the upper right panel of fig. [2] shows that 
the high-pressure ablation layer surrounding the spike is 
being compressed into a narrow region of width ~ 2/im 
and density ~ 10 g/cm 3 . The ablation also leads to signif- 
icant increase of the bubble pressure due to fluid expan- 
sion around the boundaries of the region, as can be seen 
near the spike bottom, where the expansion (and ablation) 
speed is ~ 2 - 5 x 10 5 cm/sec. Consequently, the bubbles 
on both sides of the spike are extended and deformed, and 
the fundamental spike becomes even narrower. The bub- 
bles are now nearly stationary. 

The upper right panel of fig. [2] shows that second be- 
comes strongly compressed. The secondary spikes can be 
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Fig. 2: (Color online.) Amplitudes (ablation layer displace- 
ment) of the density spikes and bubbles for the WP regime 
(a = 0.7 and /? = 0.4). For clarity the center of the density pat- 
tern has been moved to y = Ai/2. Upper left: the fundamental 
(or master) spike and bubble pattern at t — 7.78 ns. Upper 
right: the secondary density spikes (around y — 0, 12^m) at 
t — 8.10 ns. Lower left: the mass flux pu at t — 8.56 ns. Lower 
right: the fundamental spike ruptures at t — 8.73 ns. 
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Fig. 3: (Color online.) The amplitudes of the density spikes 
versus time. Upper left: the fundamental mode and the sec- 
ond, third, and forth harmonics for the WP case. Upper right: 
the same for the SP1 case. Lower panel: the lengths of the 
fundamental and second spikes, and their ratio. 




Fig. 4: (Color online.) Evolution of the fundamental spike and 
bubble for the SP1 regime (a = 5.0, /3 = 1.0). The initial 
modulation wavelength is 12/^m. Left: t = 8.93 ns, right: t = 
9.67 ns. 



attributed to the second harmonic (^(i)) having the wave- 
length Ai/2. The corresponding rate of fluid displacement 
can be larger than that of the fundamental spike and thus 
compete for mass and energy feed with it. In fact, at 
t > 8.1 ns one finds &(t) > £i(£) in the left panel of fig. [31 
i.e., the secondary spikes are overtaking the fundamental 
spike. 

As already invoked a priori in the interpretation of the 
simulation results, the evolution of the ARTI can be un- 
derstood in terms of the nonlinear interaction among the 
fundamental and the secondary spikes associated with the 
harmonics of the fundamental excitation. For convenience, 
we shall adopt the single-mode approach often used for 
considering the nonlinear RTI 2,3,16,19 25,. Generation 
of higher harmonics in the region of interest and their in- 
teractions shall be investigated. An unstable initial (the 
fundamental) single mode perturbation with wave num- 
ber fci = 2tt/Ai can excite harmonic modes with wave 
numbers k n = nki, where n = 2,3,... 2 , 3] [14] . The to- 
tal surface displacement of the interface can thus be ex- 
pressed as £(j/,i) = Y, n $n(y,t)cos(k n y). In the weakly 
nonlinear stage, £(y,t) can be represented by the first 
three harmonic modes [19] . As the instability grows, more 
modes are included. Accordingly, we first describe the 
structure of the ARTI modes by the fluid displacement 

e(y-Ai/2,t) = En £»(|/-Ai/2,*)cosMi/-Ai/2)], where 
n = 1, 2, 3, ... with n = 1 being the fundamental mode (to 
maintain symmetry in the analysis, in this representation 
the y coordinate has been displaced by Ai/2). Accord- 
ingly, one obtains £ so (±Ai/2,t) =-&(*) + 6(i) - &(*) + 
U - &(±Ai/6,t), where &(±Ai/M) - [&(*) - &(*) " 
%4,(t)]/2 — &(£) is the displacement at the bottom of the 
bubble, or£ sa (±Ai/2,i) = -3[&(t)-&(t)-$i(t)]/2. Since 
& - 6, we have U(Ai/2,t) - 3&(t)/2 > 0. 

In fig. [3] the profiles of £„(£), n = 1 - 4, obtained by nu- 
merically solving the governing mode equations are shown. 
Using the value of £4 at t = 8.1 ns from the upper- left 
panel, we find £, sa — 1.2/xm (for the same instant) from the 
lower panel. This result agrees well with the correspond- 
ing value 1.3/Ltm (for the secondary spikes in fig. [2] upper 
right panel) from the simulations. That is, by including 
four harmonics in the theory the obtained maximum dis- 
placement £ so of the secondary spikes agrees with that of 
the simulations. 

For £ sa < and £1 > £2 > •••, the distinct maximum 
in the linear ARTI growth rate becomes smoothed out, 
and there is no secondary spike excitation. Instead, at 
t = 8.1 ns new cavitation (see fig. [2] upper- right panel) 
near the secondary spikes appear. The cavitation and 
the fundamental bubble coalesce into an irregular, trans- 
versely widened, large bubble between the fundamental 
and the secondary spikes. The fundamental spike is fur- 
ther compressed and becomes still narrower. At t = 8.56 
ns, new eddies of the mass flux are generated between 
the fundamental and secondary spikes, as indicated by 
arrows in the lower-left panel of fig. [2] The mass flux 
pushes against the waist of the fundamental spike, such 
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that the density of upper region of the fundamental spike 
increases to p > 7.5 - 10 g/cm 3 and that of the lower 
part is less. The upper part is accelerated upward and 
becomes a jet. The secondary spike grows faster than the 
fundamental one. The ratio R = £ sa (±Ai/6, t)/£ ma (0, t), 
where £ mo (0,t) = ^2 n ^ n {t) ( n =1 ~ 4), is a measure of 
the displacement of the fundamental spike-bubble pattern 
at the center. It increases gradually, as shown by the red 
and green curves in the lower panel of fig. El That is, the 
secondary spikes acquire more fluid mass than the funda- 
mental spike. The plasma at the foot of the fundamental 
spike also feeds the secondary spikes, making them grow 
faster, as shown in the lower left panel of fig. [2 where the 
arrows indicate the direction and magnitude of the mass 
flux. Eventually, at t = 8.73 ns the fundamental spike 
ruptures into two parts, as observed in the simulation re- 
sult (fig. [21 lower right panel). Simulations also show that 
for the parameters considered the threshold above which 
the fundamental spike will eventually rupture is R ~ 0.20. 

If the initial modulation wavelength is larger than 
Ai=12 /im, we found by including in the analysis harmon- 
ics higher than the fourth that at least 14 higher harmon- 
ics are excited and the first 12 modes have initial growth 
rates larger than that of the fundamental mode. That 
is, many different modes can be excited simultaneously 
and the long-wavelength fundamental mode in our model 
is not preferred (even though it is started with a finite 
amplitude) . The overall grow rate is less than the preced- 
ing cases and the evolution behavior of the SP and WP 
regimes can no longer be distinguished, as can also be 
seen from the linear growth-rate versus wavelength curves 
(see fig. [U right panel). In fact, the ARTI growth rates 
for longer wavelength perturbations varies slowly with the 
wavelength. In the simulation result we see that two other 
spikes (corresponding to the additional shorter wavelength 
perturbations) appear in the original cavity region. The 
fundamental spike eventually ruptures due to lack of mass 
feed resulting from competition for the latter among the 
three spikes. 

As can also be seen in the right panel of fig. Q] the lin- 
ear growth rate in the SP regime is smaller than that of 
the WP regime, and it has no sharp maximum. This is 
because the SP regime corresponds to longer preheating 
length L and lower Atwood number. For the SP case of 
a — 5.0 and /? = 1.0 (denoted by SP1), the initial pertur- 
bation and wavelength remain the same as that without 
preheating. The cutoff wavelength is about 2.8/im and 
the linear growth rate has no sharp maximum. For the 
a = 7.0 and /3 = 1.4 case (denoted by SP2), the cut- 
off wavelength is about 3.6/im. The growth rate versus 
wavelength curve is relatively smooth and the local maxi- 
mum is not sharp. Examining the results of the SP1 case 
in more detail, we found that except in the early linear 
growth phase, no harmonic mode grows faster than the 
fundamental mode, as shown in the right panel of fig. [3] 
This can be attributed to the fact that the input energy is 
spread over several modes, instead of concentrating only in 



the mode with the highest (sharply peaked) growth rate, 
as in the WP case. In fact, in the SP regime the growth of 
the fundamental and the secondary harmonics in the non- 
linear stage is similar. As a result, the fundamental spike 
does not rupture. Instead, it evolves into a long jet that 
spurts towards the higher-temperature underdense plasma 
at an average speed of ~ 1.8 x 10 6 cm/sec, as shown in fig. 
|U Such long jets have been frequently observed in experi- 
ments (3l2l)] as well as astrophysical clouds [?3]. However, 
it should be cautioned that when several harmonics are 
excited simultaneously, the single-mode assumption can 
become invalid and other modes that are not harmonics 
of the fundamental mode can also be involved. In this case 
turbulent mixing is expected to dominate. 

In summary, we have shown that the evolution of ARTI 
in the presence of preheating can be interpreted in terms 
of the generation and interaction of harmonics of the fun- 
damental perturbation. In agreement with that from ex- 
isting numerical simulations [ISHMl 22 , it is found that 
the initial growth of the fundamental spike-bubble system 
depends on the ablation length. However, the evolution of 
the ARTI also depends significantly on the preheating pa- 
rameters. In fact, there can be two regimes of ARTI evolu- 
tion. The WP regime is characterized by a distinct peak in 
the wavelength dependence of the linear growth rate and 
rapid excitation of the corresponding, namely the second, 
harmonic. In the nonlinear stage of the evolution of the in- 
stability the fundamental and the second harmonic modes 
compete with each other. Eventually the second harmonic 
dominates and the fundamental spike ruptures. In the SP 
regime, on the other hand, the linear growth is weaker and 
not sharply peaked, so that many harmonics are excited. 
In the nonlinear evolution the fundamental spike remains 
dominant and evolves into a long jet that spurts towards 
the higher-temperature underdense blow-off region 2,3. 
Both scenarios of the ARTI evolution have often been ob- 
served in the laser-plasma compression experiments as well 
as in astrophysical plasmas ITIf5T-f51 f2"5"Il2l)] . 
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